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We numerically study the behavior of two-dimensional skyrmions in the presence of a quasi-one- 
dimensional sinusoidal substrate under the influence of externally applied dc and ac drives. In the 
overdamped limit, when both dc and ac drives are aligned in the longitudinal direction parallel to the 
direction of the substrate modulation, the velocity-force curves exhibit classic Shapiro step features 
when the frequency of the ac drive matches the washboard frequency that is dynamically generated 
by the motion of the skyrmions over the substrate, similar to previous observations in superconduct¬ 
ing vortex systems. In the case of skyrmions, the additional contribution to the skyrmion motion 
from a non-dissipative Magnus force shifts the location of the locking steps to higher dc drives, and 
we find that the skyrmions move at an angle with respect to the direction of the dc drive. For a 
longitudinal dc drive and a perpendicular or transverse ac drive, the overdamped system exhibits 
no Shapiro steps; however, when a finite Magnus force is present we find pronounced transverse 
Shapiro steps along with complex two-dimensional periodic orbits of the skyrmions in the phase- 
locked regimes. Both the longitudinal and transverse ac drives produce locking steps whose widths 
oscillate with increasing ac drive amplitude. We examine the role of collective skyrmion interactions 
and find that additional fractional locking steps occur for both longitudinal and transverse ac drives. 

At higher skyrmion densities, the system undergoes a series of dynamical order-disorder transitions, 
with the skyrmions forming a moving solid on the phase locking steps and a fluctuating dynamical 
liquid in regimes between the steps. 

PACS numbers: 75.70.Kw,75.25.-j,75.47.Np 


I. INTRODUCTION 

Phase locking or synchronization effects can arise in 
coupled oscillators when the different frequencies lock to¬ 
gether over a certain range of parameter space, an effect 
that was first reported by Huygens for the synchroniza¬ 
tion of pendulum clocks^. Phase locking has been ex¬ 
tensively studied for numerous dynamical systems rang¬ 
ing from a pair of coupled oscillators to an entire cou¬ 
pled oscillator array^i^. A single particle moving over 
a tilted one-dimensional washboard potential can also 
experience phase locking when an additional ac driving 
force is applied. The substrate periodicity produces in¬ 
trinsic periodic modulations of the particle velocity in 
the absence of an ac drive which increase in frequency 
as the magnitude of the tilt or dc drive increases. Ad¬ 
dition of an external fixed-frequency ac drive produces 
locking regimes in which the average dc velocity remains 
constant even as the magnitude of the dc drive is in¬ 
creased. The same picture can be applied to Joseph- 
son junctions, where the analog of a velocity-force curve 
is the voltage-current curve, which exhibits a series of 
phase locking regions called Shapiro steps under an ap¬ 
plied ac drive for single junctions^i^ and coupled arrays 
of junctions^. One of the hallmarks of Shapiro steps is 
that the step width oscillates as a function of the ac 
drive amplitude^^^. Shapiro step phenomena also arise 
in dc and ac driven charge density waves^^^, spin den¬ 
sity wavesiS, and Frenkel-Kontorova models consisting of 
commensurate or incommensurate arrangements of par¬ 


ticles moving over ordered or disordered substratesii^H. 
In the case of vortex motion in type-II superconductors, 
Martinoli et al. reported the first observation of Shapiro 
steps for dc and ac driven vortices interacting with a 
periodic one-dimensional (ID) substrate created by peri¬ 
odic thickness modulations of the samplei^d^, while sim¬ 
ilar effects were observed for vortices driven over 
or two-dimensional (2D)i^d^ periodic substrates. More 
recently, Shapiro steps have been found for ac and dc 
driven colloidal particles moving over a quasi-lD peri¬ 
odic substrate^S. 

Shapiro steps can also occur when a lattice of collec¬ 
tively interacting particles moves over a random substrate 
under combined dc and ac drives. Here, the effective elas¬ 
tic coupling between the particles comprising the lattice 
generates an intrinsic washboard frequency that can lock 
to the applied ac driving frequency. Such steps have been 
studied for vortices moving over random disorder— ^2^ or 
through confined channel geometries^^. For particles con¬ 
fined to 2D and moving over a quasi-lD substrate, both 
the ac and dc drives must be applied in the same direction 
to produce Shapiro steps; however, for vortices moving 
over 2D periodic or egg-carton substrates, it is possible 
to obtain what are called transverse phase locking steps 
when the ac drive is perpendicular to the direction of the 
dc drive^l^S^. These phase locking steps are distinct from 
Shapiro steps, and their widths grow quadratically with 
increasing ac amplitude rather than showing the oscilla¬ 
tory behavior associated with Shapiro steps. Phase lock¬ 
ing effects can also occur for overdamped particle motion 
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in 2D periodic systems under combinations of two per¬ 
pendicular ac drives, producing localized and delocalized 
motion as well as rectification effects^^— . 

In systems such as vortices and colloidal particles, an 
overdamped description of the equations of motion is ap¬ 
propriate. In contrast, the skyrmions that were recently 
discovered in chiral magnets have particle-like proper¬ 
ties and many similarities to superconducting vortices, 
but have the important distinction that there is a strong 
non-dissipative Magnus force in their motion^^^— . The 
skyrmions can be set into motion by an applied cur¬ 
rent and are observed to have a very small depinning 
threshold^^— in part because the effectiveness of 
the Magnus force can be up to ten times stronger than 
the dissipative force component. The Magnus force intro¬ 
duces a velocity component of the skyrmion that is per¬ 
pendicular to the direction of an imposed external force, 
so a skyrmion deflects from or spirals around an attrac¬ 
tive pinning site rather than moving directly toward the 
potential minimum as would occur for systems governed 
by overdamped dynamics^i^i^^— . Since skyrmions are 
particle-like objects, many of their dynamical proper¬ 
ties can be captured using a point particle model based 
on a modified Theile’s equation that takes into account 
repulsive skyrmion-skyrmion interactions, the Magnus 
force, damping, and substrate interactions^^*^. Such 
as approach has been shown to match well with micro- 
magnetic modeling^ of the depinning of skyrmions in 
periodic^ and random pinning arrays^. Particle-based 
skyrmion models were used to describe the motion of 
skyrmions interacting with single pinning sites^i^ as 
well as skyrmion motion in confined regions^. Since 
skyrmions can easily be driven with an applied external 
drive they potentially open a new class of experimentally 
accessible dynamical systems where the Magnus force has 
a dominant effect. It should be possible to create vari¬ 
ous types of potential energy landscapes for skyrmions 
through techniques such as thickness modulations, pe¬ 
riodic applied stain, controlled irradiation, or spatially 
periodic doping. An open question is how known phase 
locking phenomena would be affected by the inclusion of 
a Magnus force, and whether new types of phase lock¬ 
ing effects might appear that are absent in overdamped 
systems. Skyrmions also have potential for various spin- 
tronic applications^ which would require the skyrmions 
to move in a controlled manner, so an understanding of 
skyrmion phase locking dynamics could be useful for pro¬ 
ducing new methods for precision control of skyrmion 
motion. 

In this work we examine Shapiro steps for skyrmions 
moving over a quasi-lD periodic washboard substrate 
similar to the geometry considered by Martinoli et al. 
for vortices moving over substrates with a periodic thick¬ 
ness modulationi^^— . We consider the motion of single 
skyrmions and collectively interacting skyrmions over a 
periodic substrate, as illustrated in Fig. [TJ Here the lon¬ 
gitudinal or cc-direction is aligned with the direction of 
the periodicity of the substrate, while the y-direction cor- 



FIG. 1: (Color online) Skyrmions (red dots) at a density of 
Pa = 0.001 on a periodic quasi-lD substrate with Ap = 1.0. 
The darker regions are potential maxima and the lighter re¬ 
gions are potential minima, while lines indicate the skyrmion 
trajectories, (a) For an ac drive = 1.0 applied in the 
longitudinal or i-direction at amjoLd = 1.0 and = 0, 
the skyrmions oscillate in ID paths at a 45° angle to the x 
axis, (b) For an ac drive 1^“° = 0.75 applied in the transverse 
or j/-direction with = 0 at amlctd = 0.0 or the over¬ 
damped limit, the skyrmions move in ID paths along the y 
direction, (c) The same as in (b) with cxmloLd = 1.0, where 
the skyrmions form elliptical 2D counterclockwise orbits. 

responds to the transverse direction. In this geometry 
in the overdamped limit, Shapiro steps occur when the 
dc and ac (F“'^) driving forces are both applied 

along the a;-direction. When a finite Magnus force is 
present, we find that phase locking continues to occur; 
however, the net skyrmion motion is rotated at an angle 
with respect to the substrate periodicity direction. As 
the magnitude of the Magus force prefactor is increased, 
the width and number of steps for a fixed dc drive grad¬ 
ually decreases, and the steps shift to higher values of 
In the overdamped limit, if the ac drive is applied 
in the transverse or y-direction, no Shapiro steps occur, 
but for a finite Magnus force a new set of Shapiro steps 
can arise, with the width and the number of resolvable 
steps in the velocity-force curve increasing with increas¬ 
ing Magnus force prefactor. On a phase-locked step, the 
skyrmion motion takes the form of intricate periodic 2D 
orbits. As a function of the ac amplitude, we show that 
the widths of the phase locking steps for both longitu¬ 
dinal and transverse ac driving exhibit the oscillatory 
behavior associated with Shapiro steps. 

We note that the Shapiro steps we observe are differ- 
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ent than the previously studied phase locking motion for 
purely dc-driven skyrmions moving over a 2D periodic 
substrate^. In the latter case, the phase locking was 
associated with a directional locking effect in which the 
skyrmion motion locks to certain symmetry directions of 
the substrate potential. In the present study there is no 
phase locking without an ac drive. 

II. SIMULATION AND SYSTEM 

We consider a 2D system of size L x L with periodic 
boundary conditions in the x and y directions contain¬ 
ing Ns skyrmions at a density of ps = Ng/L'^. Single 
[Ng = I) or multiple skyrmions interact with a quasi-lD 
periodic sinusoidal potential with a periodicity direction 
running along the x direction, as illustrated in Fig. [T] 
The equation of motion for a single skyrmion i with ve¬ 
locity Vi = dxi / dt moving in the x — y plane is 

aw. + a^z X V. = Ff + Ff + F'"" + F“. (1) 

Here is the location of the skyrmion and ad is the 
prefactor of the damping force that aligns the skyrmion 
velocity in the direction of the net external forces. The 
second term is the Magnus force with prefactor am , which 
rotates the velocity into the direction perpendicular to 
the net external forces. In order to maintain a constant 
magnitude of the skyrmion velocity we impose the con¬ 
straint a^ + a^ = 1 and vary the relative importance of 
the Magnus force to the damping force by changing the 
ratio Om/arf. In the overdamped limit am = 0.0, while 
for skyrmions amjctd can be ten or large r'^^d'^ . 

The skyrmion-skyrmion interaction force is F®'* = 
where % = \ri - rj\, Xij = {vi - 
rj)/Rij, and Ki is the modified Bessel function. This in¬ 
teraction is repulsive and falls off exponentially for large 
Rij. For most of this work we remain in the limit where 
skyrmion-skyrmion interactions are weak so that we can 
consider the dynamics of a single skyrmion; however, we 
show that most of our results are robust under the inclu¬ 
sion of skyrmion-skyrmion interactions. The substrate 
force F®^ = —VU{xi)x arises from a washboard poten¬ 
tial 

U{x) = Uo{cos{2TTXi/a)) (2) 

where Xi = Vi ■ x, a is the periodicity of the substrate, 
and we define the substrate strength to be Ap = 2 ttUo / a. 
Unless otherwise noted, we take Ap = 1.0. The dc 
driving term is slowly increased in magni¬ 

tude to avoid any transient effects. The ac driving term 
is either F^ = cos(a;t)x for transverse driving or 
F“ = cos(a;t)y for perpendicular driving. 

We measure the time-averaged skyrmion velocities in 
the X direction (14) = ' ^)/NgUja and y- 

direction (I 4 ) = X]i=i ■ y)/NgOJa. Here, due to 

the periodicity of the substrate, phase locked steps occur 


when the skyrmions travel integer multiples of the sub¬ 
strate periodicity na during each ac drive cycle, allowing 
us to label the steps n = 0 for the pinned phase and 
n = 1,2... for the higher order steps. We focus on the 
two ac frequencies w = 8 x 10“^ inverse simulation time 
steps for the longitudinal ac driving and w = 1.6 x 10“^ 
inverse simulation time steps for the transverse ac driv¬ 
ing, and use a substrate lattice constant of a = 3.272. 

We use two different driving protocols as illustrated in 
Fig. [ 1 ] For longitudinal driving, we have 

pdrive ^ pac cos(a;t)x, (3) 

corresponding to the conditions under which Shapiro 
steps arise for an overdamped system. For transverse 
driving, we have 

pdrive ^ pdc^ p pac cos(cdt)y, (4) 

which would produce no Shapiro steps in the overdamped 
limit. 

In Fig. [Tfa) we show the skyrmion trajectories for 
= 1.0, = 0.0, amloLd = 1.0, and a skyrmion den¬ 

sity of Ps = 0.001. In this case the skyrmions are pinned 
and form a triangular lattice that is commensurate with 
the substrate. The ac drive causes the skyrmions to os¬ 
cillate in the potential minima; however, their motion 
is not strictly in the ^-direction but is tilted at an an¬ 
gle oi 9 = 45° with respect to the x direction due to 
the Magnus force, which induces a velocity component 
perpendicular to the ac driving direction. In the ab¬ 
sence of a substrate, a dc or ac drive applied in the 
x-direction causes the skyrmions to move at an angle 
9 = arctan(am/ad) with respect to the driving direction, 
so that in the overdamped limit of am = 0.0 the skyrmion 
moves parallel to the direction of the net external driv¬ 
ing force. In Fig. [Jb), we rotate the direction of the ac 
drive to be in the transverse direction with F“° = 0.75 
and F'^^ = 0 for a sample with amjotd = 0.0. In this case 
the skyrmion motion follows strictly ID paths aligned 
with the j/-direction that pass through the potential min¬ 
ima of the substrate. For amjctd = 1-0, as shown in 
Fig. He), the skyrmions rotate in counterclockwise ellip¬ 
tical patterns, showing that the Magnus force can induce 
x-direction motion even when the drive is applied only 
in the y-direction. In the absence of the substrate the 
ac drive would produce only ID trajectories at an angle 
with respect to the y-axis. This highlights the fact that 
the Magnus force affects how the skyrmions move when 
interacting with forces induced by the substrate. 


III. LONGITUDINAL AC DRIVING 

We first consider the case illustrated in Fig. Ha) of ac 
driving in the longitudinal direction. We conduct a series 
of simulations for increasing amjOLd and focus on the sin¬ 
gle skyrmion limit. In general we find that the Shapiro 


4 



FIG. 2: (14) (upper blue curves) and (14) (lower red curves) 

vs F'^‘^ for the system in Fig.[T](a) in the single skyrmion limit 
at = 1.0. (a) In the overdamped limit of amlctd = 0, 
( 14 ) = 0 and a series of steps appear in (14) indicating phase 
locking, (b) At am/ctd = 0.58, (I 4 ) is finite, (c) am/ad = 
3.042 and (d) am/o-d = 9.962 show the increase of skyrmion 
motion in the direction transverse to the substrate and the 
shift in the locking phases. 


steps we observe remain robust when finite skyrmion- 
skyrmion interactions are included; however, additional 
features can arise for varied fillings when the skyrmion 
structure is incommensurate with the substrate, as we 
discuss in Section V. In Fig. ^a) we plot (14) and (Vy) 
versus for the system in Fig. UKa) at = 1.0 in the 
overdamped limit of arajo-d = 0. Here (V),) = 0 while 
( 14 ) shows a series of steps indicative of the phase lock¬ 
ing. These features are similar to those observed for other 
overdamped systems moving over quasi-lD periodic sub¬ 
strates such as vortices in type-II superconductors mov¬ 
ing over quasi-lD substrate modulations. In Fig. HKb), 
when arajo-d = 0.58, both (Vy) and (14) are finite and 
have a ratio of |(14)/(14)| ~ 0.58. Here the phase lock¬ 
ing is still occurring, but the intervals of in which 
the phase locking steps appear are shifted. Figure HKc) 
shows that at ctmlctd = 3.042, both |(I4)| and some of 
the step widths have increased in size, and there are no 
clear regions between the steps where no phase locking is 
occurring. In Fig. [2Kd), at am/oid = 9.962, there is only 
a single phase locking step. 

To more clearly demonstrate the behavior of the steps 
for varied Om/ad, in Fig. [3Ka) we plot (14) versus F‘^‘^ 
for amjoid ranging from 0.0 to 11.147, with the evolution 
of the first three locking steps n = 1 to 3 highlighted. 



FIG. 3: (a) (14) vs F^^ at Ap = 1.0 for am/o-d = 0.0 

(brown), 0.577 (light blue), 0.98 (dark purple), 1.33 (light 
purple), 2.06 (dark orange), 3.042 (light orange), 4.0 (dark 
red), 4.92 (light red), 7.0 (dark green), 8.407 (light green), 
9.962 (dark blue), and 11.147 (black), from left to right. Here 
( 14 ) exhibits quantized values corresponding to specific steps, 
(b) The corresponding values of (Fy) vs which contains 
steps that are not quantized. 



FIG. 4: The regions of phase locking for the n = 0 to n = 8 
steps as a function of and am/od. The width of the steps 
is reduced and the steps shift to higher values of F'^° with 
increasing Qm/od. 


For a given value of n, the step in (Vr) has a fixed value 
regardless of the choice of Om/ckd, and each step shifts 
to higher values of F'^^ with increasing amlo-d- The cor¬ 
responding (Vy) versus plot in Fig. [3Kb) shows that 
the steps in (I 4 ) are not quantized in integer multiples 
of 2tt jau. The quantization of the (14) arises from the 
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FIG. 5: (14) vs F'*" for a^/a^ = 9.962. (a) = 2.4. (b) 

= 4.2. (c) Fr = 6.2. 


periodicity of the substrate in the a;-direction, and since 
the y-direction has no periodicity, there is no quantiza¬ 
tion of (Vy). In Fig. |4] we highlight the evolution of the 
widths of the n = 0 through n = 8 steps as a function of 
and amjoLd at = 1.0. At Umlctd = 0, the largest 
number of phase locking steps can be resolved. We ob¬ 
serve two general trends as am/o-d increases. First, for 
n > 3, the widths of the locking regions decrease and 
the intervals of over which the locking occurs shift 
to higher values of with the magnitude of this shift 
increasing with increasing n. Second, the width of the 
n = 1, 2, and 3 steps initially increases for increasing 
dm!Cid before reaching a maximum and then decreasing 
again. The width of the n = 0 step reaches a maximum 
with increasing am!ctd and then saturates. The shift in 
the locations of the phase locking regions arises because 
the angle at which the skyrmions move with respect to 
the x-axis increases with increasing amlcx-d, causing the 
skyrmions to spend larger intervals of time interacting 
with the repulsive portion of the substrate potential. As 
a result, higher values of must be applied to cause the 
skyrmion to translate in the x-direction at larger am!ctd- 
We next determine if the phase locking steps at high 
Magnus force prefactor are of the Shapiro type. In Fig.[S] 
we plot {Vx) versus F‘^‘^ for amjctd = 9.962 at F®"^ = 2.4, 
4.2, and 6.2 to show the variation in the widths of the n = 
0, n = 1, and n = 2 steps. In Fig. [5] we plot the widths 
Ao and Ai of the n = 0 and n = 1 steps, respectively, 
versus F“'^. Each step shows the characteristic oscillation 
expected for Shapiro steps, where the width of step n 
is proportional to \Jn{F^‘^)\, where is the nth-order 
Bessel function^. The solid lines in Fig. inKa,b) are fits to 
I Jo I and |Ji|, respectively. The higher order steps obey 



FIG. 6: (a) The width Aq of the n = 0 step vs for the 
system in Fig. [S] at am/oLd = 9.962. The solid line is a fit to 
the I Jo I Bessel function, (b) The width Ai of the n = 1 step 
vs Fx’^ for the same system. The solid line is a fit to the |Ji| 
Bessel function. In each case the width of step n shows an 
oscillation of the form of the Bessel function | Jn(F“'^)|, which 
is characteristic of Shapiro step phase locking. 



FIG. 7: (Vx) vs F'^'^ for a system with dc driving in the 
x-direction and ac driving Fy^ = 1.0 in the y-direction. (a) 
At OLmloLd = 0.0, there are no steps in (W). (b) At amlo-d = 
0.436, steps are present, (c) am/ad = 2.06. (d) amjad = 
3.04. (e) arnlccd = 4.0. (f) amlcid = 9.962. 


similar fits. This indicates that in the Magnus-dominated 
limit, Shapiro step phase locking is occurring. 


IV. TRANSVERSE AC DRIVING 

We next consider the case illustrated in Fig. [TKb,c), 
where the ac drive is applied transverse to the direction 
of the substrate potential. In the overdamped limit of 
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o-mlo-d = 0, such a drive causes the skyrmion to oscil¬ 
late in the y-direction as shown in Fig. [TJb), and when a 
finite dc drive is applied in the longitudinal direction, a 
single washboard oscillation frequency in the ^-direction 
is generated by the motion of the skyrmion over the pe¬ 
riodic substrate. Since only one frequency is present, 
there is no coupling between two frequencies, so mode 
locking does not occur. When the Magnus force is fi¬ 
nite, the transverse ac drive induces an oscillating ve¬ 
locity component in the longitudinal or x-direction as 
well as in the y-direction, as illustrated in Fig. me), so 
that it is possible for the dc-induced washboard frequency 
to couple to the transverse ac frequency and generate 
a transverse Shapiro step. In Fig. [3 we plot (14) vs 
for a single skyrmion moving with = 1.0. At 
o-rajo-d = 0, shown in Fig. m^)) there are no steps in 
( 14 ), indicating the lack of phase locking, while the cor¬ 
responding (Fy) = 0. Depinning occurs at the threshold 
value Fc of = Ap = 1.0. Figure [TDb) shows that at 
oimlotd = 0.436, the depinning threshold has dropped 
substantially to Fc = 0.4 and a series of steps are now 
visible for 0.4 < F'^'^ < 2.0, indicating that phase lock¬ 
ing is occurring. For amlc^d > 0.0, (Vy) is finite and the 
|(Fy)| versus F'^^ curve has exactly the same form as (14) 
versus F'^^, but the magnitude of |( 14 )| is multiplied by 
drnlctd- In Fig. m^jd) we plot (F,) versus F'^'^ for sam¬ 
ples with amjaa = 2.06 and 3.04, respectively. Here, 
the widths of the locking steps increase with increasing 
o-mlo-d and the step locations are shifted to higher val¬ 
ues of F'^^. In samples with amlctd = 4.0 and 9.962, as 
shown in Fig.[7Ke,d), respectively, the steps extend out to 
larger values of F'^^, and the non-phase locking regions 
between the steps are also extended. The steps in (14) 
once again occur at quantized values of naa;/27r due to 
the periodicity in the x-direction, while the steps in (Fy) 
do not have quantized values. 

In Fig. [5] we plot the location of the upper edge of 
the n = 0 step as a function of F^^ and amlctd for the 
system shown in Fig. [S] with Ap = 1.0. This is equivalent 
to the threshold depinning force Fc. Here F'^^ jAp = 1.0 
at amjcLd = 0.0 and it decreases to zero at amlctd = 
1.226. There is a local maximum in F^'^lAp at amlctd = 
1.55, followed by another minimum near amlad = 2.2 
and a broad plateau for higher values of amlad- This 
oscillatory behavior in the n = 0 step width is absent 
for longitudinal ac driving, as shown in Fig. [5] where Fc 
exhibits only monotonic behavior. The dips and maxima 
in Fc for the transverse ac driving are associated with 
transitions in the shape of the skyrmion orbits during a 
single ac drive cycle for increasing amlad- 

In Fig. |9] we illustrate the skyrmion trajectories in 
a subsection of the system on the n = 0 step at the 
points labeled (a-d) in Fig. [5] Figure EDa) shows that 
for amlad = 0 and F‘^‘^ = 0.1, the skyrmion moves in 
a ID path in the y-direction along the potential mini¬ 
mum. At amlad = 0.75 and F"^^ = 0.2, in Fig. IHDb), the 
skyrmion forms an elliptical orbit that is confined within 
a single potential trough. On the local maximum in the 



FIG. 8: The location of the upper edge of the n = 0 step as 
a function of and otmloLd for the system in Fig. [5] with a 
transverse ac drive of F^'^ = 1.0. Here there are several local 
minima and maxima that are associated with changes in the 
skyrmion orbits, as shown in Fig. [9] at the points marked a-d. 



FIG. 9: Skyrmions (red dots), potential maxima (darker re¬ 
gions), potential minima (lighter regions), and skyrmion tra¬ 
jectories (lines) in a portion of the system in Fig. [8] along 
the n = 0 step at the points labeled (a-d) in Fig. [8] (a) 

At a-mlctd = 0.0 for F'^'^ = 0.1, there is ID motion in 
the y-direction. (b) amlctd = 0.75 at F'^'^ = 0.2. (c) At 
amlad = 1.54 and = 0.1, the skyrmion moves between 
two potential minima, (d) At amlad = 4.92 and F'*'^ = 0.2, 
the skyrmion moves between three potential minima. 
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FIG. 10: (a) Evolution of the regions in which the n = 0, 1, 2, 
3, 4, 5, 6, 7, 8, 9, 10, and 11 steps (from bottom to top) appear 
as a function of and otmloid for Fy‘^ = 1.0. Increasing the 
Magnus force produces enhanced phase locking, (b) A blowup 
of panel (a) in the region of small amjad showing that the 
steps vanish as amjad goes to zero. 


n = 0 step marked point c in Fig. |S1 at arajoLd = 1.54 
and = 0.1, the skyrmion forms a more complicated 
2D orbit that has three lobes. In a single ac drive cycle 
the skyrmion translates back and forth by two substrate 
lattice constants. The dip in Fc at am/ad = 1.226 shown 
in Fig. [5] corresponds to the point at which the skyrmion 
orbit transitions from being confined in one potential 
minimum to traversing two potential minima. Above 
the second local minimum at amloLd = 2.2 in Fig. |S1 
the skyrmion orbit becomes even more complex, as illus¬ 
trated in Fig.[3](d) for amloid = 4.92 and = 0.2. The 
skyrmion now moves between three substrate potential 
minima in a single ac drive cycle. The local minimum in 
the n = 0 step width at amictd = 2.2 then corresponds 
to the transition in the skyrmion motion from travers¬ 
ing two substrate minima to traversing three substrate 
minima. For higher values of ttmlcid, additional minima 
in Fc could occur that would be correlated with orbits 
traversing four or more substrate minima. We expect 
that additional substrate minima would be resolvable in 
samples with a smaller substrate lattice constant a. 

In Fig. [Tol a) we highlight the regions of phase locking 
as a function of F‘^'^ and amjad for steps n = 0 through 
n = 11 for the system in Fig. 0 When amloLd = 0, all 
the steps with n > 1 vanish, as illustrated in Fig. fTOlb) 
where we plot the regime 0 < amloid < 1.0. As the Mag¬ 
nus force increases, a larger number of steps can be re¬ 
solved. In general, the step widths increase with increas¬ 
ing arn/od! however, certain steps such as n = 1, 2, and 
3 show step width oscillations. In the case of longitudinal 
ac driving, the skyrmion orbits along the different lock¬ 





FIG. 11: Skyrmions (red dots), potential maxima (darker re¬ 
gions), potential minima (lighter regions), and skyrmion tra¬ 
jectories (lines) for the system in Fig. (9] (a) n = 1 orbit at 
Om/od = 0.75 and F'^'^ = 0.6. (b) n = 2 orbit at amloid = 
0.75 and = 0.7. (c) n = 2 orbit at amjad = 1.55 and 
Fdc _ Q n = 1 orbit at amjad = 2.06 and = 0.3. 

Here the skyrmion translates by 2a in the positive x direction 
followed by a in the negative x direction for a net transport 
by a distance a in the i-direction during each ac cycle, (e) 
n = 3 orbit at amjad = 2.06 and F'^‘^ — 0.67. (f) n = 1 orbit 
at amjad = 5.92 and F^'^ = 0.67. 


ing steps are always ID in nature. In contrast, the orbits 
are much more complicated for transverse ac driving. In 
Fig . fTTl a) we show the n = 1 skyrmion orbit from Fig. [TUI 
at amjad = 0.75 and F'^'^ = 0.6. The skyrmion translates 
in the positive cc-direction and negative y-direction, mak¬ 
ing an angle close to 9 — arctan(am/od) = 36.9° with the 
x-axis. During a single orbit the skyrmion passes through 
a loop and translates by one lattice constant in the x- 
direction. Figure [TTT b) illustrates the n = 2 orbit for 
amjad = 0.75 at F'^‘^ = 0.7, where the skyrmion trans¬ 
lates two lattice constants in the x-direction per ac cycle. 
On the n = 2 step for amjad = 1.55, and = 0.5, 
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FIG. 12: (a) (14) vs for = 1.0, = 9.962, 

and Ap = 0.5 (black), 2.0 (green), 4.0 (blue), and 7.0 (red), 
(b) The evolution of the n = 0, 1, 2, and 3 step widths as a 
function of F'^’^ and Ap for the system in panel (a). 


shown in Fig. HHc), the skyrmion moves at a steeper an¬ 
gle oi 9 = 57.1° from the x-axis. In Fig. fTTT d). which 
shows the n = 1 orbit at amjctd = 2.06 and = 0.3, 
during a single ac drive cycle the skyrmion initially moves 
2a in the positive x-direction before moving a in the neg¬ 
ative x-direction, producing a net translation in the x- 
direction of a distance a per ac cycle. Figure fTlT e) shows 
the amloid = 2.06 system in the n = 3 orbit at = 67, 
where the skyrmion translates by 3a in a single ac cycle. 
On the n = 1 step at amldd = 4.92 and F'^^ = 0.67, 
plotted in Fig. ITIT f). the skyrmion moves 3a in the pos¬ 
itive x-direction during the first portion of the ac drive 
cycle followed by 2a in the negative x-direction during 
the second portion of the ac drive cycle, producing a net 
translation of a in the x direction during a single ac cy¬ 
cle. We observe similar orbits for the other values of n, 
and find that the net angle of the skyrmion motion with 
respect to the x axis increases with increasing amjad- 



FIG. 13: (a) The width Aq of the n = 0 step vs F^° for the 
system in Fig. [T2]at Om/Qd = 9.962 and Ap = 1.0. The solid 
line is a fit to the |Jo| Bessel function, (b) The width Ai of 
the n = 1 step vs F“° for the same system. The solid line is 
a fit to the |Ji| Bessel function. 

the locking regions oscillates with increasing Ap, and for 
Ap > 8.0 all the locking phases shift linearly to higher 
values of F'^'^ with increasing Ap. The step width oscil¬ 
lations arise due to variations in the number of potential 
minima through which the skyrmion orbit passes during 
a single ac drive cycle, similar to what was observed for 
fixed Ap and varied amjad- This result shows that the 
transverse phase locking is a generic feature that appears 
in both the strong and weak substrate regimes, and that 
it is more pronounced for stronger substrates. 

We next examine the dependence of the step widths 
at a hxed Ap on the ac driving amplitudes, as shown in 
Fig. [Ql where we plot Aq and Ai versus F“° for Ap = 1.0 
and amjctd = 9.962. The solid lines are fits to A„ oc 
|J„(F“°)|, indicating that the transverse phase locking 
steps are also of the Shapiro step type, similar to the 
longitudinal phase locking steps. 

V. COLLECTIVE EFFECTS 


A. Dependence on Substrate Strength and ac 
Amplitude 

We next consider the effect of the substrate strength 
on the transverse locking steps at amlctd = 9.962 and 
F“° = 1.0. In Fig. [UKa) we plot (14) versus F'^'^ for 
Ap = 0.5, 2.0, 4.0, and 7.0. At the lower values of Ap, 
the phase locking steps decrease in width, and the steps 
completely vanish when Ap = 0. This is highlighted in 
Fig. HHb) where we plot the widths of the n = 0, 1, 2, 
and 3 steps as a function of and Ap. The width of 


We next consider assemblies of interacting skyrmions 
for the system shown in Fig. [TJ In general, when the 
skyrmion density is commensurate with the substrate and 
the skyrmions can form a triangular lattice, skyrmion- 
skyrmion interactions cancel and we find the same types 
of phase locking observed in the single skyrmion sys¬ 
tems. For incommensurate fillings where dislocations 
are present or when the skyrmion structure becomes dis¬ 
torted or anisotropic in the pinned phase, we find that it 
is possible for additional fractional phase locking to occur 
between the integer phase locking steps. These fractional 
locking steps occur when a portion of the skyrmions are 
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FIG. 14: (14) vs F‘^‘^ at OLmlotd = 2.0. (a) ac driving in the 
a;-direction with = 1.0 for a single skyrmion (dark blue 
line) and a sample containing multiple skyrmions at a density 
of ps = 0.04 (light orange line), showing that fractional phase 
locking steps can arise, (b) The same for ac driving in the 
y-direct ion at Fy’^ = 1.0. 


locked to step n and the remainder of the skyrmions are 
locked to step n — 1. In Fig. IHT a) we plot (14) versus 
for a system with ac driving in the cc-direction at 
c^mlctd = 2.06 and = 1.0 to compare the results for 
a single skyrmion with a system at a skyrmion density 
of Ps = 0.04. There are no fractional steps in the single 
skyrmion system; however, when interacting skyrmions 
are present we find fractional steps njm, where n and 
m are integers. Figure lldf b) shows the same system for 
ac driving in the y-direction, where the same types of 
fractional steps arise. The fractional steps appear at 
incommensurate fields when it is possible to have two 
effective particle species in the system. One species is 
commensurate and the other is associated with intersti¬ 
tials, dislocations, or vacancies. In overdamped systems 
such as superconducting vortices moving over 2D peri¬ 
odic substrates, similar integer steps for individual or 
non-interacting vortices appear at commensurate match¬ 
ing fillings while additional fractional locking steps arise 
at non-matching fieldsi^. 

At much higher skyrmion densities and for sufficiently 
strong substrate strengths, the pinned skyrmion struc¬ 
tures become highly anisotropic due to the confinement 
in the ID pinning rows. In the moving phase just above 
depinning, the effectiveness of the pinning is partially 
reduced and the repulsive skyrmion-skyrmion interac¬ 
tions favor a more uniform structure. The competition 
between skyrmion-skyrmion and skyrmion-substrate in¬ 
teractions produces a series of order-disorder transitions 
in the moving state. On the phase-locked steps, the 



FIG. 15: (a) (K) vs F'^’^ at = 1.0 and = 2.06 

for a sample with a skyrmion density of pa = 0.4. (b) The 
fraction of six-fold coordinated particles Pe vs for the 
same system showing that along the phase-locked steps the 
skyrmions form a much more ordered state. 


skyrmions form an ordered moving anisotropic lattice 
and travel in a synchronized fashion, while between the 
phase locking steps the skyrmions adopt a more isotropic 
or liquid like configuration. 

In Fig. fTKI a) we plot (14) versus F‘^‘^ for a sample with 
= 1.0, amjctd = 2.06, and a skyrmion density of 
Ps = 0.4, showing the n = 1, 3, and 4 phase locking steps. 
Figure fTSf b) illustrates the corresponding fraction of six¬ 
fold coordinated skyrmions Pq = N~^ ~ Zi), 

where Zi is the coordination number of skyrmion i ob¬ 
tained from a Voronoi construction. On the phase lock¬ 
ing steps, Pq increases to Pq = 0.92, while between the 
steps Pq ~ 0.5 on average and shows strong fluctuations. 
In Fig. fTBI a.b) we show the real space locations of the 
skyrmions and the corresponding structure factor S{k) 
on the n = 1 step at F'^‘^ = 0.25 from Fig. [151 The 
skyrmions are all moving together and form a partially 
ordered but anisotropic lattice. Even though the system 
is anisotropic, most of the skyrmions have six neighbors, 
so that Pq Ri 0.9. Figure [TBI c.d) shows the same sample 
at = 0.6, corresponding to the non-phase locking re¬ 
gion labeled b in Fig. [T51 Here the skyrmions form a dis¬ 
ordered structure that is less anisotropic than the phase 
locked state. We observe similar sets of dynamical order- 
disorder transitions between step and non-step regions 
for increasing and find similar effects for ac driving 
in the a;-direction. Studies in overdamped systems of col¬ 
lections of interacting vortices also show that the vortices 
are more ordered and exhibit suppressed noise fluctua¬ 
tions in a phase locked region22i2^. At higher the 
effectiveness of the substrate gradually diminishes, the 
phase locking steps disappear, and the skyrmions can re- 
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FIG. 16: (a,c,e) The real space positions of the skyrmions 

from Fig. 1151 and (b,d,f) the corresponding structure factors 
S'(fc). (a,b) The n = 1 phase locked step at F'^'^ = 0.25 from 
the point labeled a in Fig. IlSl ai shows a partially ordered 
anisotropic structure. (c,d) On the non-step region at = 
0.6 labeled b in Fig. I15l ai. the skyrmions form a disordered 
liquid like structure. (d,e) On a non-step region at — 9.0, 
the system forms a moving lattice. 


These results show that Shapiro steps for skyrmions in¬ 
teracting with a periodic substrate are a robust feature 
that occurs for a variety of skyrmion densities and sub¬ 
strate strengths. The change in the skyrmion lattice 
structure as the system passes in and out of phase locked 
states as a driving current is swept could be observed 
using neutron scattering or noise measurements. 

VI. SUMMARY 


We have analyzed Shapiro steps for skyrmions inter¬ 
acting with periodic quasi-one-dimensional substrates in 
the presence of combined dc and ac drives, with a specific 
focus on the role of the Magnus force in the dynamics. 
When the dc and ac drives are both applied in the lon¬ 
gitudinal direction, which is aligned with the substrate 
periodicity, phase locking occurs, and as the role of the 
Magnus force increases, the phase locking steps gradu¬ 
ally reduce in width and shift to higher values of the 
driving force. The skyrmions move at an angle to the 
direction of the external dc drive that increases as the 
contribution of the Magnus force increases. When the 
ac drive is applied perpendicular to the dc drive and the 
substrate periodicity direction, there is no phase lock¬ 
ing in the overdamped limit; however, when there is a 
finite Magnus force, phase locking can occur. On the 
phase locked steps the skyrmions move in intricate two- 
dimensional periodic orbits. We map out the evolution of 
the phase locked regions for the transverse and longitudi¬ 
nal ac driving for varied contribution of the Magnus force, 
ac driving amplitudes, and substrate strength. When 
collective interactions between skyrmions are introduced, 
fractional Shapiro steps can appear. For strong substrate 
strengths and higher skyrmion densities, both longitu¬ 
dinal and transverse phase locking steps occur that are 
associated with dynamically induced transitions between 
an ordered anisotropic solid on the steps to a fluctuat¬ 
ing liquid state in the non-phase locked regimes. Such 
transitions could be observed with neutron scattering. 


order into a more uniform moving crystal state as shown 
in Fig. lTSl' e.fl at F'^'^ = 9.0. Similar dynamical reordering 
to a triangular lattice for high drives has been observed 
for skyrmions interacting with random pinning^ as well 
as for vortices driven over random pinning arrays^^^. 
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